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Abstract 

A statistical model of self-organization in a generic class of one-dimensional nonlinear Schrodinger (NLS) 
equations on a bounded interval is developed. The main prediction of this model is that the statistically pre- 
ferred state for such equations consists of a deterministic coherent structure coupled with fine-scale, random 
fluctuations, or radiation. The model is derived from equilibrium statistical mechanics by using a mean-field 
approximation of the conserved Hamiltonian and particle number {L^ norm squared) for finite-dimensional 
spectral truncations of the NLS dynamics. The continuum limits of these approximated statistical equilib- 
rium ensembles on finite-dimensional phase spaces are analyzed, holding the energy and particle number at 
fixed, finite values. The analysis shows that the coherent structure minimizes total energy for a given value 
of particle number and hence is a solution to the NLS ground state equation, and that the remaining energy 
resides in Gaussian fluctuations equipartitioned over wavenumbers. Some results of direct numerical inte- 
gration of the NLS equation are included to validate empirically these properties of the most probable states 
for the statistical model. Moreover, a theoretical justification of the mean-field approximation is given, in 
which the approximate ensembles are shown to concentrate on the associated microcanonical ensemble in 
the continuum limit. 

PACS: 05.20.-y; 05.45.-a; 52.35.Mw 
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1 Introduction 

The appearance of macroscopic organized states, or coherent structures, in the midst of turbulent small-scale 
fluctuations is a common feature of many fluid and plasma systems [Q. Perhaps the most familiar example 
is the formation of large-scale vortex structures in a turbulent, large Reynolds' number, two-dimensional fluid 
^. A similar phenomenon occurs in slightly dissipative magnetofluids in two and three dimensions, where 
coherent structures emerge in the form of magnetic islands with flow |^] . In the present work, we shall be 
concerned with another nonlinear partial differential equation whose solutions in the long-time limit tend to 
form large-scale coherent structures while simultaneously exhibiting intricate fluctuations on very fine spatial 
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scales. Namely, we shall consider the class of one-dimensional wave systems governed by a nonlinear Schrodinger 
equation of the form 

#t+^..+/(iV'nv^ = o, (1) 

for various nonlincaritics / chosen so that the dynamics is nonintegrable and free of wave collapse. 

Numerical simulations of the NLS equation (]l|) in a bounded spatial interval with either periodic or Dirichlet 
boundary conditions demonstrate that, after a sufficiently long time, the field ip evolves into a state consisting of 
a time-periodic coherent structure on the large spatial scales coupled with turbulent fluctuations, or radiation, 
on the small scales 0, For a focusing nonlinearity / (/(a) > 0, /'(a) > for a > 0) the coherent structure 
takes the form of a spatially localized, solitary wave. This organization into a soliton-like state is observed 
from generic initial conditions for the focusing NLS equation At intermediate times a typical solution 
consists of a collection of the soliton-like structures, which as time evolves undergo a succession of collisions, or 
interactions. In these collisions the larger of the solitons increases in amplitude and smaller waves are radiated. 
This interaction of the solitons continues until eventually a single soliton of large amplitude survives in a sea of 
small-scale turbulent waves. (Even though (|l|) does not possess solitons in the strict sense unless = IV'Pj 

it does have spatially localized solitary wave solutions in the generic focusing case. As a matter of convenience, 
and in keeping with the terminology in |§, 0], we therefore simply refer to the coherent states that emerge in 
the simulations in this case as solitons). Figure (1) below illustrates this behavior for the particular saturated 
focusing nonlinearity /(IV'P) = I^P/ (1 + IV'P)- For a defocusing /, coherent structures also emerge and persist, 
but they are less conspicuous since they are not spatially localized. The reader is referred to for a detailed 
description of numerical investigations of the long-time behavior of the system (Q) for various nonlinearities /. 

Theoretical arguments and numerical experiments suggest that these coherent structures for the NLS 
dynamics correspond to ground states of the system. That is, a coherent structure assumes the form 
ipix, t) = exp(— iAi), where 4>{x) is a solution of the ground state equation 

(t>.. + fm^)<p+x<p^o, (2) 

In turn, these ground states (j) are the minimizers of the Hamiltonian 

HW = l f l^.l^dx-^ I F{\t^\^)dx (3) 

subject to the constraint that the particle number 

^W^ljj^l'dx (4) 

be equal to given initial value iV°. We shall refer to the first and second terms in as the kinetic energy 
and the potential energy, respectively. The potential F is defined by F{a) = f{a')da' . The parameter A 
arising in (||) is the Lagrange multiplier in the variational principle: minimize i?(<^) subject to N((j)) = N^. 
For focusing nonlinearities /, the ground states (j) are spatially localized. By contrast, for periodic boundary 
conditions and for defocusing nonlinearities /, the ground states are spatially uniform; this can be proved by 
a straightforward application of Jenson's inequality when the potential F is strictly concave. Nevertheless, 
independent of the properties of the ground states, the coherent structures are expected to be constrained 
minimizers of the Hamiltonian given the particle number, whenever these minimizers exist. 

Yankov and collaborators |[ 0] have suggested that the tendency of the solution of the NLS system in 
the focusing case to approach a coherent soliton state coupled with small-scale radiation can be understood 
within a statistical mechanics framework. They argue that the solutions to the ground state equation that 
realize the minimum of the Hamiltonian subject to the constraint on the particle number are "statistical 
attractors" to which the solutions of the NLS equation (|^) tend to relax. In this argument, the process which 
increases the amplitude of the solitons as the number of solitons decreases is thought to be "thermodynamically 
advantageous", in the sense that it increases the "entropy" of the system. This entropy seems to be directly 
related to the amount of kinetic energy contained in the radiation, or the small-scale fluctuations of the field 
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ip, so that the transference of kinetic energy to the fluctuations is accompanied by an increase in the entropy. 
Similar ideas have been expounded by Pomeau ||l0|| who has argued, based on weak turbulence theory, that 
the defocusing cubic NLS equation = "IV'P) in a bounded two-dimensional spatial domain should 

exhibit the tendency to approach a long-time state consisting of solution of the ground state equation plus 
small-scale radiation. 

Our primary purpose in the present paper is to construct a statistical equilibrium model of the coherent 
structures and the turbulent fluctuations inherent in the long-time behavior of system . With this model we 
seek to translate the various intuitive ideas about the balance between order and disorder in the NLS dynamics 
into explicit, verifiable calculations. In particular, we provide strong support to the notion set forth in ^, 0, |l^ 
that the ground states, which minimize the energy for a given particle number, are statistically preferred states. 
In addition, we furnish a definite meaning to the concept that random fluctuations absorb the remainder of the 
energy. 

We derive our equilibrium statistical model from a mean-field approximation of the conserved Hamiltonian 
H and particle number N for a finite-dimensional truncation of the NLS equation (^. This approximation 
relies on the fact that the fluctuations of tJj about its mean, the coherent structure, become asymptotically 
small in the continuum limit, provided that the mean energy and mean particle number remain finite as the 
number of modes in the spectral truncations goes to infinity. In constructing the mean-field theory, therefore, 
we choose that probability measure on the finite-dimensional phase space which maximizes entropy subject to 
approximated constraints on the mean energy and particle number. Specifically, we approximate the means of 
the particle number and the potential energy term in the Hamiltonian in terms of the mean state ("0). 

This approach gives rise to a Gaussian model, which is accordingly easy to analyze. Thus, we find that to 
any initial value N*^ of the particle number there corresponds a coherent structure, which is precisely the mean 
state, or most probable state, for the statistical model. Moreover, we show that this mean state is indeed a 
solution to the ground state equation that minimizes the Hamiltonian over all states with the same particle 
number iV°. Furthermore, we clarify the role of fine-scale fluctuations by demonstrating that the gradient field 
ipx has finite variance at every point. Thus, we see that the difference between the initial energy H*^ and the 
energy of the mean state, H{{il})), resides entirely in the kinetic energy of infinitesimally fine-scale fiuctuations. 
This result, which makes precise the ideas of Yankov et al. |^] and Pomeau |ic|] , explains how it is possible 
for an ergodic solution to the Hamiltonian system to approach a solution to the ground state equation 
without violating the conservation of energy or particle number. 

In the present paper the principal justification for our mean-field approximation is the evidence of extensive 
numerical simulations, which show that the fluctuations in ^ go to zero over longer and longer time intervals 
1^, ^, Ij. Nevertheless, we also include an a posteriori justification of this approximation by proving that, in 
the continuum limit, the probability measures defining the mean-field theory with approximately conserved 
quantities H and N concentrate about the phase space manifold on which H — and N = . In other 
words, we establish a form of the equivalence of ensembles, showing that our mean-field ensemble becomes 
equivalent in a certain sense to the microcanonical ensemble, which is a measure concentrated on the manifold 
H = H'^,N = iV°. Besides connecting our mean-field approach to the standard principles of equilibrium 
statistical mechanics pT| , this result shows that the mean-field theory is asymptotically exact in the continuum 
limit. This property of the mean-field ensembles is a common feature of turbulent continuum systems in which 
the mean values of conserved quantities are fixed while the number of degrees of freedom goes to infinity [p^ , [l^ . 

These topics are taken up after a quick review of NLS theory and previous work on invariant measures for 
NLS. 



2 The NLS Equation and its Conserved Quantities 

The NLS equation (0) describes the slowly- varying envelope of a wave train in a dispersive conservative system 
and so arises in many branches of physics. It models, for example, gravity waves on deep water |p^ , Langmuir 
waves in plasmas [^|, pulse propagation along optical fibers [p^ , and self- induced motion of vortex filaments 
p^ , [l^ . In particular, the cubic NLS equation, which corresponds to the nonlinearity /(|V-'P) — ilV'P: has 
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garnered much attention in the mathematics and physics Hterature. On the whole real line, or on a periodic 
interval, the cubic NLS equation is completely integrable via the inverse scattering transform ^ . The NLS 
equation with any other nonlinearity or boundary conditions, however, is not known to be integrable. 

In the present article we are interested in the NLS equation (|l|) as a generic model of nonlinear wave 
turbulence. For this reason, we wish to consider nonlinearities / and boundary conditions such that the dynamics 
is well-posed for all time and yet is nonintegrable. For instance, a natural class of prototypes of this kind consists 
of those NLS equations with bounded nonlinearities / on a bounded spatial interval with Dirichlet boundary 
conditions. Bounded nonlinearities arising in physical applications are often called saturated nonlinearities, of 
which — IV'P/Cl + IV'P) and /(|V'P) = 1 ~ exp(— are examples. These amount to corrections to the 

focusing cubic nonlinearity for large wave amplitudes, and they have been proposed as a models of nonlinear 
self- focusing of laser beams , propagation of nonlinear optical waves through dielectric and metallic layered 
structures , and self- focusing of cylindrical light beams in plasmas due to the ponderomotive force ||2^ . Such 
nonlinearities certainly satisfy the requirements set forth by Zhidkov to guarantee well-posedness of the NLS 
equation (|l|) in the space for Dirichlet boundary conditions [Q. 

In light of these properties, we shall adopt the NLS equation with homogeneous Dirichlet boundary conditions 
on a bounded interval fl = [0, L] and with a (focusing or defocusing) bounded nonlinearity / as the generic 
system throughout our main development. Precisely, we impose on / is the following boundedness condition 

sup(|/(a)| + |(l + a)/'(a)|) < ex.. (5) 

In Section 9 we indicate how our main results can be extended to a natural class of unbounded nonlineari- 
ties. Also, when it is convenient for comparison with some numerical simulations, we can instead use periodic 
boundary conditions. The necessary modifications, which are straightforward, are left to the reader. 
The NLS equation d) may be cast in the Hamiltonian form 

i , 6H 

where -0* is the complex conjugate of the field tp and H is the Hamiltonian defined by (||). For our class of 
generic NLS systems, the total energy H and the particle number N , or norm squared, are the only known 
dynamical invariants for (^). These invariants play the leading role in the statistical equilibrium theory and the 
values and , derived from a given initial state, say, alone determine the statistical ensemble. 

For later use, we write the NLS equation (|l|) as the following coupled system of partial differential equations 
for the real and imaginary components u{x,t) and v{x,t) of the complex field ip{x,t): 

ut + v^^ + fiu'^ + v'^)v = 0, vt-u^x- fiu'^ + v'^)u = 0. (6) 

Similarly, the Hamiltonian can be written in terms of the fields u and v as 

H{u, v) ^ K{u, v) + e(u, v) (7) 

with kinetic and potential energy terms 

K{u, w) = ^ I vl) dx , Q{u, f ) = / F{u'^ + w^) dx , 



and the particle number is given by 



N{u, w) = i / {u^+ v^) dx . (8) 



In the same way, the ground state equation (g) for (/) = u -|- iv can be expressed as the following coupled system 
of differential equations: 

Uxx + f{u'^ + v'^)u + Xu = 0, Vxx + f{u'^ + v'^)v + Xv = . (9) 
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3 Invariant Measures for NLS 



Several authors have constructed invariant measures for the NLS equation (yj). Zhidkov [g4| considered nonlin- 
earities / satisfying the boundedness condition (||) and showed that the Gibbs measure 

Pp{dudv) = exp {-'!3H{u, v)) \Y du(x)dv{x) , (10) 

for fixed /3 > is normahzable and has a rigorous interpretation as an invariant measure for the NLS system on 
the phase space L^. Similar results have been obtained by Bidegary These results hold for both periodic 
and Dirichlet boundary conditions. However, Pfj as constructed by Zhidkov fails to account for the invariance 
of the particle number (^) under the NLS dynamics. 

Earlier, Lebowitz et al. j26| considered the problem of constructing invariant measures for NLS on the peri- 
odic interval with the focusing power law nonlinearities — ■ For such nonlinearities the Hamiltonian 
H is not bounded below, and therefore the measure Pp defined by ( |lO| ) is not normahzable. As an alternative, 
Lebowitz et al. conditioned the Gibbs measure on the particle number and studied the measures 

P^,Aro {dudv) = Z-^ exp{-f3H{u, v))I{N{u, v) < N°) Y[ du{x)dv{x) , (11) 

xen 

where I{A) is the characteristic function of the set A. They showed that P^jyo is normahzable for 1 < s < 4 
for any iV'^, and for s = 4 if is sufficiently small. For s > 4, though, Pfj^^a fails to be normahzable. 

The normalizability of the conditioned Gibbs ensemble (|l]) is closely related to the well-posedness of the 
NLS equation as an initial value problem . Indeed, for the focusing power law nonlinearities it can be shown 
that the NLS equation (Ej) with periodic or Dirichlet boundary conditions is well-posed, or free of wave collapse, 
in the Sobolev space H (ft) as long as s < 4, while blow-up can occur from smooth initial conditions when 
s > 4 [^6[ . Bourgain proved that under the conditions on s and -/V° considered by Lebowitz et al. 

the measures Pp. no are rigorously invariant under the NLS dynamics and he thereby obtained a new global 
existence theorem for the initial value problem with these nonlinearities. Analogous results have also been 
established by McKean 

A common feature of all these invariant measures for the NLS equation is that for any fixed (3 > the 
kinetic energy (1/2) /(u^ + u^) dx is infinite with probability 1, and hence the Hamiltonian is also infinite with 
probability 1 ||2^. This difficulty occurs whenever classical Gibbs statistical mechanics is applied to a system 
with infinitely many degrees of freedom and is related to the well-known Jeans ultraviolet catastrophe lll|| . 
Consequently, these ensembles are appropriate when the typical state of the continuum system is such that 
each of its modes contains a finite energy. They are not appropriate, on the other hand, to situations in which 
a typical state of the system is realized by an ergodic evolution from an ensemble of initial conditions having 
finite energy and particle number iV'', both of which partition amongst the many modes of the system. 

Since the purpose of the current paper is to construct ensembles which provide meaningful predictions about 
the long-time behavior of solutions of the NLS equation (Q), we are interested in statistical equilibrium ensembles 
which, in the continuum limit, have finite mean energy and mean particle number. In view of this goal, we 
scale the inverse temperature parameter /3 in the Gibbs measure with the number of degrees of freedom n and 
thereby maintain the mean energy at a finite value as n — > oo. The continuum limit we obtain under this scaling 
is degenerate in certain ways, but this is hardly surprising given that the known invariant measures for NLS 
are supported on fields with infinite kinetic energy. On the other hand, the analysis of the scaled continuum 
limit is greatly simplified by the fact that the random fluctuations of (u, v) at each point tend to zero. This 
fortuitous property implies that a natural mean-fleld approximation is asymptotically exact. Moreover, the 
mean-field ensemble concentrates on the microcanonical ensemble for the invariants H and N in the scaled 
continuum limit. For this reason, we can avoid the canonical ensemble constructed from H and N, which fails 
to be normahzable for focusing nonlinearities because the linear combination H — XN always has a direction in 
which it goes to —oo, as we shall demonstrate in Section 5. 
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4 Discretization of NLS 

We now introduce a finite-dimensional approximation of the NLS equation (|l|) with homogeneous Dirichlet 
boundary conditions on the interval — [0,L]. For this purpose we shall use a standard spectral truncation, 
even though many other approximation schemes could be invoked. Let ek{x) = \/'2/L sm(y/Xi^.x) and = 
(fc7r/L)^,fc — 1,2,..., be the eigenfunctions and eigenvalues of the operator — on with homogeneous 
Dirichlet boundary conditions. The eigenfunctions form an orthonormal basis for the real space L'^{fl). For 
given n, let X„ = span {ei, . . . , e„} C L^, define u^") and w*^"^ in X„ by 

n n 

u^''\x,t) = J2Mt)ekix), v^''\x,t) = J2Mt)ekix), (12) 

k=l fe=l 

and set ip'^^H^i ^) = u^"'>{x, t) + iv'^"'\x, t). We consider the following evolution equation for t/j^"): 

z^^") + vi".^ + p(") (/(iV'^"' nv^^")) - , (13) 

where P*^"^ denotes the orthogonal projection from L^(r2) onto X„. Equation ( p^ is clearly a spectral truncation 
of the NLS equation (^ . It is equivalent to the following coupled system of ordinary differential equations for 
the Fourier coefhcients Uk and Vk,k = 1, . . . ,n: 

^ - XkVk + ^ fiiu^'-^r + {v^-'^f)v^''^eu dx = 0, (14) 

^ + XkUk ^ /((«("))' + (v^-^fW-^ek dx = 0, (15) 

where we have used the fact that {ek)xx — —Xkek- 

Whether viewed as the evolution of (w("\i;^"^) on X„ x X„ or in terms of the Fourier coefficients {uk,Vk) 
on i?^", it can be shown ||2^, ^ that this finite-dimensional dynamical system has Hamiltonian structure, with 
Hamiltonian 7J„ = Kn + Qn, where 



if„(u("),z;(")) = I [jui-yf + {vi-yf]dx ^ if^ Afe(^.^ (ig) 
is the kinetic energy, and 
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e„ (..("), !;(")) = / P((^^("))' + (t^("))')dx, (17) 



is the potential energy. The functionals i?„, if„ and 0„ are just the restrictions to Xn x X„ of the functionals 
H,K and 0, which are defined by (0). The Hamiltonian Hn is, of course, an invariant of the dynamics. The 
particle number 



2 7o 2 

•'^ fc=i 



is also conserved by this dynamics, as may be verified by direct calculation. Also, 7V„ is the restriction of the 
functional N, which is defined by (Sh, to the space Xn x X„. 

The Hamiltonian system (p^)-(15|) as a dynamical system on R^" satisfies the Liouville property In 
other words, the Lebesgue measure 11^=1 dukdvk on R^" is invariant under the phase flow for this dynamics. 

5 Mean-Field Ensembles for NLS 

We now proceed to construct a statistical model that describes the long-time behavior of the truncated NLS 
system (|l^)-([l5|). According to the fundamental principles of equilibrium statistical mechanics, the Liouville 
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property and the ergodicity of the dynamics provide the usual starting point for such a description, and the 
microcanonical ensemble is the appropriate statistical distribution (that is, probability measure on phase space) 
with which to calculate averages for an isolated Hamiltonian system |o). In the specific system under 
consideration, the natural canonical random variables are Uk and Vk, which comprise the phase space i?^". 
The random fields ^("^(a;) = '^Ukek{x), = "^Vkekix), and = u'^^\x) + as well as 

the functions iJ„ — H {u^'^\ v^'^^ ) and iV„ — 7V(m("\ w^")), are determined by these canonical variables. The 
microcanonical ensemble for the truncated NLS system therefore takes the form 

n 

Puo^Noidudv) = W-^5{Hr,- H°)5{Nn- N^)Wdukdvk, (19) 

fc=i 

where W — W{H^ , N^) is the structure function or normalizing factor. Under the ergodic hypothesis, expecta- 
tions with respect to this ensemble equal long-time averages over the spectrally-truncated dynamics from initial 
conditions with prescribed, finite values and N'^ of the invariants and Nn- 

While the microcanonical ensemble is a well-defined invariant measure for the spectrally-truncated dynamics, 
it is a cumbersome to analyze, and the calculation of ensemble averages for finite n is not feasible. In statistical 
mechanics the usual procedure to overcome this difficulty is to introduce another invariant measure, the canonical 



ensemble, and prove that in the limit as n — )■ oo, it becomes equivalent to the microcanonical ensemble 11 
For the NLS system the canonical Gibbs ensemble is 



Pl3,x{dudv) = Z-^exp{~p[Hn~ XNn]) Y[dukdvk, (20) 



fc=i 



where /3 and A are chosen such that the averages of Hn and iV„ with respect to Pp^x equal H° and iV°, re- 
spectively, and Z = Z(f3, A) is the partition function or normalizing factor. Unfortunately, for general focusing 
nonlinearities /, even those satisfying the boundedness condition (|^), this canonical measure is not normal- 
izable. To see this degeneracy, let us consider a nonlinearity / that is both positive and strictly increasing 
and an associated ground state (j) that is a nontrivial solution to (^. Then, for any positive scale factor cr, a 
straightforward calculation yields the identity 

(H-XN){acf>) = Y / - M^^dx -\[ F{a^\cl,\-') dx 



dx r [/(a)-/(|0n]da 
SI Jo 



It follows from this expression that, as a increases to infinity, the value of _ff — XN at acf) tends to — oo. The 
same reasoning applies to the discretized system. Consequently, the function iJ„ — AiV„ in the Gibbs weight has 
a direction in which is goes to —oo, resulting in the divergence of the partition function Z in (po|). In essence, 
this divergence reflects the fact that such a ground state is a critical point, but not a minimizcr, of H — XN. 
Locally, the first variation SH — X6N vanishes at (p, but the second variation 6^H — X6^N is positive only on 
those variations 5(p for which SN = 0. When (50 = (1 — cr)0 for ct near 1, and thus 6N ^ 0, the corresponding 
second variation of H — XN is negative. The above calculation shows that this local behavior extends to a global 
degeneracy for a wide class of nonlinearities. As Lebowitz et al. discuss in this defect of the canonical 
ensemble is even more severe when the nonlinearity / is unbounded. 

These considerations compel us to find another way to approximate the microcanonical ensemble and thus 
to derive a tractable statistical equilibrium model. First, we note that since the microcanonical measure (^) 
depends on u*^"^ and u^"-* through i7„ and N^ only, it is invariant with respect to multiplicative constants of 
the form e'^. In view of this phase invariance, we can factor the random field in the form e*^0^"\ where 9 
is uniformly distributed on [0, 27r] and is independent of (/)("\ In what follows, we shall write ^^"^ = w*^"^ -|- iz;'^") 
with the phase normalization 

arg / 0(")(x)dx = 0, 
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and consider the microcanonical distribution conditioned by this constraint. It is this conditional distribution 
that we seek to approximate with a mean-field approximation. In constructing the approximation, we can 
therefore describe distributions for 0^"' that need not possess invariance under a change of phase, with the 
understanding that we can recover the random field V'^"'' by setting t/'^"-' = e*^(/)'^"^. 

The key to the mean-field construction is supplied by direct numerical simulations of the evolving microstates 
■0 governed by (|]) 0, D . Since the particle number and the Hamiltonian are well conserved in these simulations, 
we may think of these numerical experiments as realizations of the microcanonical ensemble. The simulations 
clearly demonstrate that the amplitude of fluctuations of the random field (the phase normalized field as- 
sociated with ip) at each point x G Q become small in the long-time limit, and that they contribute little to 
the conserved norm of </) ^ (also, see Figure 1 below). Moreover, this effect becomes more apparent 
as the spatial resolution of the numerical simulations is improved |^. In particular, these fluctuations become 
negligible compared to the magnitude of the coherent structure that emerges in the mean field (p, which in the 
focusing cases takes the form of a single, localized soliton. On the basis of these properties of the simulated 
microstates, we shall adopt the hypothesis that, in the limit of infinite resolution, the fluctuations of (f> are 
infinitesimal. Precisely, we shall make the following vanishing of fluctuations hypothesis: 



With respect to the phase-invariance conditioned microcanonical ensemble on 2n-dimensional phase-space, there 
holds 

lim / [Var(u(")(2;))+Var(w("'(2;))]rfa; = lim V[Var(wfe) + Var(i;fe)] = , (21) 

n — ^oo n — >oo ^ — ^ 



k=l 



where Var denotes the variance of the indicated random variable with respect to that ensemble. The first 
equality follows from the identity ^[Var(M(") (a;)) + Var(i;(")(a;))] dx = ELi[Var(Mfe) + Var(wfe)]. 

Our strategy is to use the vanishing of fiuctuations hypothesis to build a mean-field ensemble that is equiv- 



alent in an appropriate sense to the microcanonical ensemble (19) in the continumn limit as n — > cx3, but is 
easy to analyze for finite n. Rather than attempt to verify this hypothesis a priori, we shall show a posteriori 
that the resulting mean-field ensemble concentrates on the microcanonical manifold at fixed and iV°. The 
analysis of the fundamental microcanonical ensemble itself will be taken up in a subsequent publication. 

We note however that the numerical simulations of the NLS system which support the hypothesis ( pT| ) 
also strongly indicate that the fluctuations of the gradients u^^^ , wi"'' are not negligible in the limit of inflnite 
resolution. Figure 2 illustrates the typical long-time behavior of ji/'ajp for the particular nonlinearity /(|'0p) = 
|V'P/(1 + in which flnite fluctuations in the gradient tjjx persist. This property of the microcanonical 

ensemble is compatible with the vanishing of fluctuations hypothesis, as can be seen from the straightforward 
identity 

n 

[Var(4") {x)) + Var(4") {x))]dx = V [Var(Mfe) + Var(i;fe)] . (22) 
^ k=i 

In fact, as we shall see, this quantity remains finite in the continuum limit according to our statistical theory, 
and it represents the portion of kinetic energy absorbed by infinitesimally fine-scale fiuctuations. 

We now proceed to construct a mean-field ensemble for each n based on approximations derived from the 
hypothesis (pT|). This ensemble is determined by a probability density p^^^ on i?^" with respect to Lebesgue 
measure Ofc^i dukdvk- (Unlike the microcanonical ensembles, the mean-field ensembles are taken to be abso- 
lutely continuous with respect to Lebesgue measure for each n.) In order to define that p^"-* which governs the 
mean-field theory, we appeal to standard statistical-mechanical and information-theoretic principles [^H [ll| , |30|] 
and choose p*-"^ so that it maximizes the Gibbs-Boltzmann entropy functional 

^ n 

S{p) = p{ui, . .. ,Un,vi,. . . ,w„)logp(ui, . . . ,Un,vi, . . . ,Vn)Y\_ dukdvk , (23) 

k=l 

subject to some appropriate constraints. The entropy S has the well-known interpretation as a measure of the 
number of microstates {u,v) corresponding to the macrostate p ||ll|. The form of the entropy as a functional of 
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p is dictated by the Liouville property of the dynamics which requires that the entropy be relative 

to the invariant measure Ofe=i dukdvk- Alternatively, from the information theoretic point of view, maximizing 
the entropy S amounts to finding the least biased distribution on 2n-dimensional phase space compatible with 
the constraints and the uniform prior distribution Jlfc=i dukdv^ 

The constraints imposed in the Maximum Entropy Principle (MEP) are the crucial ingredient in the deter- 
mination of the ensemble p'-"-' . As is well-known, the microcanonical ensemble is produced when they are taken 
to be the exact constraints 7J„ — H'^ and Nn = N^; the canonical ensemble, which however is ill-defined in 
the focusing case, corresponds to average constraints {Hn) = H° and {Nn) = Throughout our discussion 
of (MEP) the angle brackets (•) denote expectation with respect to the admissible density p in (MEP). For 
our mean-field theory we seek constraints that are intermediate between those giving the microcanonical and 
the canonical ensembles. To this end we invoke the vanishing of fluctuations hypothesis (|2^) and impose the 
following constraints: 

1 " 

2E((^'=)' + (^'=)') = ^°- (24) 

k=l 

^EAfc((^^) + (t;^))-i f F{{u^-)f + {v^-^f)dx^H\ (25) 
k=i ''^ 

which we shall refer to as the mean-field constraints for obvious reasons. We note that these constraints involve 
only the first and second moments of the flclds u^"-* and with respect to p. 

We can motivate this choice of constraints in (MEP) by approximating the values of the conserved quantities 
Hn and 7V„ in terms of means with respect to the microcanonical ensemble. An immediate implication of ( pT| ) 



applied to the definition of Nn given in (18) is that 



N° ^ I f ((M(")(a;))2 + (t;(")(a;))2)d2; + i / [Var(M(")(a;)) + Var(z;(") (x))] 

2 .la 2 



= 7V„((u(") ),(!;("))) + 0(1) (26) 

as n — > oo. When we drop the error term in this approximation at finite n and we replace expectation with 
respect to the microcanonical distribution by expectation with respect to an admissible density p, we obtain 



the mean-field constraint (|24| ) 



Similarly, we can obtain the mean-field constraint (|25| ) from the vanishing of fluctuations hypothesis (21) by 



analyzing Hn with respect to the microcanonical distribution for large n. The kinetic energy (16) in H'^ = {Hn, 
is retained exactly as {Kn)- On the other hand, the potential energy (^^ is approximated by expanding F 
about the mean ((«(")),(«("))), which yields 

e„(w("),f(")) = -\l F((u("))2 + (t,("))2)ds 



/((w("))2 + (w("V) ((w("))(u(") - (u("))) + (w("))(v(") - («(")))) dx 
- i („H _ (^(«))) j An ,v )(^ (^(„)_^^(„)^) j rfa;, (27) 

4uvf'{u^ + v^) 2f{u^ + v^)+4v^f'{u^ + v^) J 

is the matrix of second partial derivatives of F{u^ + u^) with respect to u and v and (t2'^"\u(")) lies between 
(«("),«(")) and ((u(")), («("))). Because we are considering nonlinearities / satisfying (|^) each element of the 
matrix J(t2^"'', u*^")) is bounded over the domain Q independently of n. Therefore, taking the expectation on 



where 
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both sides of equation and noting that the second term on the right hand side of this equation has zero 
mean, we obtain 

(e„) = -i + dx + Rn, 

where the remainder i?„ satisfies 

\Rn\ < C (^[(«^"^ - ("^"^))' + («^"^ - (^^^"^»'] dx'j 

= C I [Var(u(")(x)) + Var(u(")(x))]da; 



for some constant C independent of n. The vanishing of fluctuations hypothesis ( pl| ) thus implies that i?„ 
as n — > oo, and hence that 

H° = ^EA,((^^) + («^))-i [ + + 0(1). (29) 

k=l •'^^ 

The mean-field constraint ( p5| ) results from neglecting the error term in this expression for finite n and replacing 
microcanonical expectations by expectations with respect to p. 

In summary, the mean-field theory is defined by a probability density p*^"-* on i?^" that maximizes the entropy 
S given in ( p^ ) subject to the constraints (^4|) and (|2^). These mean- field ensembles, which solve the governing 
(MEP), are analyzed in detail in the subsequent sections. Here we merely note that they enjoy some properties 
not shared by either the microcanonical or the canonical ensembles. First, unlike the microcanonical ensemble, 
the mean-field ensembles are analytically tractable because the densities p*^"' are Gaussian. This desirable 
property is a simple consequence of the fact that the constraints ( p4[ ) and (|5|) involve only the first and second 
moments of /a*-"-*. Second, in contrast to the canonical ensemble, the mean-field ensemble exists in both the 
focusing and defocusing cases. This crucial property depends on the fact that fluctuations in Nn are suppressed 
in the mean-field ensemble. 

However, because of the approximations made in deriving the constraints (E3) and (p5|), the mean-field 
ensemble p*-"-* is not an invariant measure for the truncated NLS dynamics (|l^-(]l5|) at finite n, even with the 
random phase e'^ included. Nevertheless, it becomes consistent with the fundamental microcanonical ensemble 
in the limit as n ^ oo, in the sense that (iV„) iV°, {Hn) -ff", and the variances of N„ and Hn converge to 
. These properties, which are proved in Section 8, imply that the ensembles p^"-' concentrate about the phase 
space manifold iJ„ = H'^,Nn = iV*^, and thereby justify the approximations made in deriving the mean- field 
theory. 



6 Equilibrium States for the Mean-Field Theory 

We now proceed to calculate and analyze the solutions p'"^ to the maximum entropy principle (MEP). First, 
we calculate the density p*^"^ which maximizes entropy subject to the mean-field constraints (p^-(p5|) by the 
method of Lagrange multipliers. If we denote the left-hand sides of (04) and (pq) by Nn and Hn, then we have 



SS ^ ti^"'>SNn+ (3^"''SHn, (30) 
where S denotes variation with respect to the density variable p and p^"^ and Z?^"^ are the Lagrange multipliers 



to enforce the constraints (24) and ( p5|) on Nn and _ff„. The variation of S is 

. n 

SS = - log{p)6p Y[ dukdvk , (31) 



k=i 
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for variations 5p satisfying jj^^n ^pY\1=i dukdvk = 0. Similarly, the variations of -/V„ and Hn are 



<57V„ 



n » n 

/ {{uk)uk + {vk)vk)5pY]_dukdvk, (32) 

fe=i"'-R^" k=i 



and 



(5/f„ = - / Afc (ufe + Wfe) (5/9 TT dukdvk 
k=i 2 fe=i 



using the relation F' = f. Equating the coefficients of 6p in ( ^o|) and recalling that u'"^ = X]fe=i ""/cSfe and 
^(n) _ ^^^j^tij,efe, we discover after some algebraic manipulations that the entropy-maximizing density p'") 
has the form 



p'"^(ui, . . . ,Un,Vl, . . . ,Vn) ^ Yl Pk{uk,Vk) , (34) 



where, for k = 1, . 



with 



Pk{uk,Vk) = ^ exp [{uk - {uk)f + {vk - {vk)f] I , (35) 



(«fe> = ^ /((«("H^))' + (^^(")(^)>')(^^"n^)>efe(a:) do: - Jn^^^^^) ' (36) 

(«fe> = ^ jj{{u^'^\x))^ + (z;(«)(x))2)(«(«)(a;))efc(a;)dx- . (37) 

We see immediately that mi, . . . , u„, wi, . . . , are mutually independent Gaussian random variables with 
means satisfying (|3^)-(|37|), and variances 

Var(ufe) = Var(«fc) = . (38) 

The multiplier is therefore necessarily positive. The equations (p6|)-(|3^) for the means {uk) and {vk) can 
be written as an equivalent complex equation for the mean field (0'^"') = (u^"'^) + 

+ (/(l(0(")>P)(0("^)) + A(")(0(")) = 0, (39) 

where we introduce the real parameter A*^"-* = — /i*^"'//?*^"^. Recalling that p(") is the projection onto the span 
Xn of the first n eigenfunctions ei, . . . , e„, we recognize this equation as the spectral truncation of the ground 
state equation (||) for the continuous NLS system (|]). In other words, the mean field {4>^"^) at each finite n is 
a critical point of the functional Hn — X^^^'Nn = 0, by virtue of (|3|). Thus, we draw the important conclusion 
that the mean field predicted by the statistical equilibrium theory is a solution of the ground state equation for 
the NLS system. 

Since the maximum entropy distribution p*-"-* satisfies the mean-field Hamiltonian constraint (|2^), a direct 
calculation reveals that 

^° = iEAfc(Var^.fe+Varz;fc) + if^A,((ufe)2 + («fe)2)-i /■ + 

fe=l k=l ''^ 

= ■^+i/„((u(")), («("))), (40) 
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where (38) is invoked to obtain the second equahty. In this expression the term n//?'-"'' in represents the 
contribution of fluctuations to the energy, while iJ„((u''"-'), (?;*•"-')) is the energy of the mean state, or coherent 
structure. We can rearrange (|o|) to get the following expression for /^^"^ in terms of the number of modes n 
and the energy of the mean state: 

' i/" -i/„((w(")), («(«))) ■ ^ ' 

Now we come to the following central result, which enables us to characterize completely the mean-field 
statistical ensembles p*^"^ . 

Theorem 1. The mean field pair ((u*^"^) , (u^"))) corresponding to any solution p'") of the constrained variational 
principle (MEP) is an absolute minimizer of the Hamiltonian iJ„ over all (u("\w(")) eXnX which satisfy 
the particle number constraint A'"„(u'"\ u*^"^) — N'^. 

Proof: To prove this assertion, we calculate the entropy of a solution p^"^ of (MEP). Referring to equations 
(p4)-(p3), which define any solution of (MEP) we immediately obtain 



5(p(")) = - / p("'logp(")f[dufedz;fc 

k=i 

n ^ 

= Pk{uk,Vk)\ogPk{uk,Vk)dukdvk 



E log 



2tt 

C{n) + nlog 



- 1 



fe=i 

£2 



where C(n) = n — X]fc=i log(fc27r/2) depends only on the number of Fourier modes n. Here, we have used the 
identity Afe = (kn/L)^. From this calculation and equation (^l|), it follows that 

5(/-')^C(„)+„log("''^°-^""°"""-'"»lV (42) 



Evidently, the entropy S{p^"-'^) is maximized if and only if the mean field pair ((m^")), (w'^"^)) corresponding to p*-"^ 
minimizes the Hamiltonian 7J„ over all fields (u'^"\ «(")) g X„ x X„ that satisfy the constraint Nn{u'^^\ w'"^) = 
N'^ . That such niinimizers exist is proved in the Appendix. 

The equation (p^ ) has the interesting interpretation that, up to additive and multiplicative constants, the 
entropy of the mean-field equilibrium p^"^ is the logarithm of the kinetic energy contained in the fluctuations 
about the mean state, which is the coherent structure. Because the potential energy is determined entirely by 
the mean, the difference between the prescribed total energy H'^ and the energy of the mean must be accounted 
for entirely by the contribution of the fluctuations to the kinetic energy. The same conclusion also follows from 
equation (|4C|). This result provides a new and definite interpretation to the notion set forth by Yankov et al. 
1^, |6[ and Pomeau that the "entropy" of the NLS system is directly related to the amount of kinetic energy 
contained in the small-scale fluctuations, or radiation, and that the dynamical tendency of solutions to approach 
a ground state that minimizes energy at a conserved particle number is "thermodynaniically advantageous." 

The parameter A*-"-* is determined as a Lagrange multiplier for the mean-field variational problem: minimize 
Hn subject to the constraint that iV„ = . The multiplier associated with any minimizer = (u*-"-*) -|- 

i(y^^)^ is unique, even though the minimizer itself may be nonunique. Indeed, the family of solutions e^^ {(j)^^^) 
for arbitrary (constant) phase shifts 6 exists for every such minimizer. Moreover, the ground state equation 
is a nonlinear eigenvalue equation, whose solutions can bifurcate. These bifurcations play the role of phase 



12 



transitions in the mean-field statistical theory. Some analysis of the solutions of the mean-field equation is 
provided in the Appendix. 

Once a solution pair A(")) is found, with i7„((</>("))) = i7* being the minimum value of Hn allowed 

by the particle number constraint Nn = , the "inverse temperature" /S*^"^ is uniquely determined by 

n 



= TTTT^ • (43) 



Indeed, this is a restatement of Since i/* approaches a finite limit as n oo, we find that /?*^"-' scales 

linearly with the number of modes n. This asymptotic scaling of inverse temperature, /?^"^ ~ nj3* , where (3* is 
finite in the continuum limit, is what distinguishes our mean- field theory from the statistical equilibrium theories 
discussed in Section 3. In those theories the inverse temperature is fixed and the mean energy is allowed to go 
to infinity in the continuum limit. Our rescaling of inverse temperature with n is necessary in order that the 
expectation of Hamiltonian with respect to p^"-* remain finite in the limit as n — s- cx). 

The implications of this scaling of /J^"^ are most evident in the the particle number and kinetic energy 
spectral densities. Substituting the expression ( |4^ ) for /3*^"^ into (38), we obtain the following formulas for the 
variances of the Fourier components and Ufei 

Var {uk) = Var K) = ' ^" . (44) 

From these equilibrium expressions we obtain a sharp form of the vanishing of fluctuations hypothesis (|2l|), 
which we adopted to derive the mean-fleld ensembles; namely, 

^ V[Var(iife) + Var(«,.)] = ^ = o{-) asn-xx., (45) 

using the fact that A^^ — 0{k^^) in the estimate. Thus, we have the following prediction for the particle 
number spectral density 

Uul + vl) = ^{{ukf + {v,f) + . (46) 

Since the mean-field is a smooth solution of the ground state equation, its spectrum decays rapidly in k, and 
we have the approximation ^{uf, + vj,) « {H^ — H*)/{nXk) for k >> 1/L. In the same way, the kinetic energy 
spectral density is 

^^ki{ul + vD) = l\,i{u,r + («,)') + , (47) 

2 2 n 

This expression shows that, as may be anticipated from a statistical equilibrium ensemble, the contribution to 

the kinetic energy from the fluctuations is equipartitioned among the n spectral modes. These predictions are 

tested against the results of direct numerical simulations in , where a good agreement with the mean-field 

theory is documented. 



7 Continuum Limit of Random Fields 

In this section we investigate the limits as n ^ oo of the random fields and their gradients. The 

random fields u^"-* and u^"^ will be seen to converge in a uniform way to deterministic limits u* and v* , while the 
gradients u* and v* converge in a weaker sense to random fields with finite variance but no spatial coherence. 
A simple form of the convergence of u*-"-* can be seen from the following easy estimate: 

n 

Var(u(")(a;)) = ^ Var(wfc) (efc(a;))2 
fc=i 




= O(n-i), (48) 
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and similarly for (x). Thus, the variance at each point goes to zero as n ^ oo. 

A more detailed analysis is complicated by the fact that the mean fields are known to converge 

only on appropriately selected subsequences. This result is the content of Theorem Al proved in the Appendix. 
It states that such subsequences of mean fields converge in the C^-norm to ground states {u* , v*) that solve the 
continuum equation (H); these ground states minimize the Hamiltonian H over all {u,v) G 11^(^1) that satisfy 
the particle number constraint N{u,v) — N'^. The possible nonuniqueness of solutions to the ground state 
equation necessitates the introduction of subsequences, since in general there may exist a set of limits {u*,v*) 
for a specified constraint value iV*^. Accordingly, we shall assume throughout this section that n goes to infinity 
along a subsequence for which the limiting ground state {u*,v*) exists, without relabeling n. 

The detailed properties of the continuum limit are expressed in the following theorem. 

Theorem 2. The (sub)sequence of random fields ) whose distribution is the mean-field ensemble 

with density p*-"' converges in distribution to a non-random field {u* , v* ) which minimizes the Hamiltonian 
H given the particle number constraint N — N'^. In addition, the gradients {u'i^\vi"'^) converge in the sense 
of finite-dimensional distributions to a Gaussian random field {U'{x),V'(x)) on 17 = [0, L] with the following 
properties: U' and V are statistically independent, the mean of (U',V') is {u*,v*), and U' and V each have 
the covariance function 

rrO _ rr* f 0, if X ^ y 

T*{x,y) = I 1, if a; = y and < X < L (49) 

I 2, if a; = y and x ^ 0, L 

where H* is the minimum value of H given N = N'^ . 

Proof. Recall that for each n, the mean ((u(")), (?;("))) under is a solution to ( p9[ ) which minimizes 
over Xn X X„ subject to the constraint iV„ = N^. By Theorem Al, the (sub)sequence ((u'"^), (u*^"-')) converges 
in C^(fl) to {u*,v*), where {u*,v*) is a solution of (^ that minimizes H subject to the constraint N{u, v) = N^. 

Recalling (^ and the independence of the Fourier coefficients Uk , we may use the following representation 

of 



k—l fc=l 



{uk) H , ^ , Zk 



efc 



r 

= (u(")) + ^==5]-Z,e,, (50) 

V M k=l 

where Zi , ^2 , . . . are independent Gaussian random variables with mean and variance 1 . The sum in the last 
line is the partial sum of the Karhunen-Loeve expansion of the Brownian bridge on which is known to 
converge to its limit B(x),x € f2, uniformly with probability one as n — s- oo ([^^, Theorem 3.3.3). That is, for 
each function h : C(il) — s- R which is continuous in the supremum norm, we have, with probability one, 

h I ^ —ZkBk I h{B) as n oo. (51) 

It is clear from this demonstration that for such h we have 

/i(m(")) ^ h{u*) as n ^ oo (52) 

with probability one, since along this (sub)sequence the means converge uniformly and Z?^"-* goes to infinity. 
Restricting to bounded h and taking averages, we have 

(/i(u("))) h{u*) as n ^ oo, (53) 

which is the definition of convergence in distribution of the random fields u*-"-* to u* . The argument that w'^") 
converges to v* is identical. 
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The analysis of the gradient fields Ux^^ and is simpler. First note that they are independent and differ 

(n) 

only in their means, and consequently it is sufficient to consider Ux , say. From the expression 

n 

fc=l 

it is evident that ui"-* is a continuous Gaussian random field whose mean («!"■') converges uniformly to u* as 
71 ^ oo. Its covariance is then calculated to be 

1 " 1 

k—1 

^ [D^inix + v)/L) + Dn{^{x - y)/L) - 2] , (55) 



where 

Dn{e) = 



2L/3(") 

( sin(n+l)9 



for e ^ 0,±27r,±47r, . 



2n+l, for 61 0, ±27r, ±47r, . 



is the Dirichlet kernel The asymptotic behavior /J*^"' '--^ n/{H'^ — H*) implied by ( p3| ) combined with the 
properties of Dn{9) then produce the desired result that, as n ^ cxd, the covariance function Tn{x, y) converges 
pointwise to the function r*(a;, y) given in the statement of the theorem. 

Now let U' be the Gaussian random field with mean u* and covariance T* . For convergence of finite- 
dimensional distributions, it remains to show that for each m and any . . . , Xm in the joint distribution on 
i?™ of the vector (^^"^(xi), . . . , ui^\xm)) converges weakly to the joint distribution of (C/'(xi), . . . , U'{xm)) as 
n — > oo. But this follows by the Levy Continuity Theorem |Q, since the characteristic function of the former 
vector converges to the characteristic function of the latter vector, simply because the means and covariances 
converge. The argument that wi"^ converges to V' is similar, as is the extension to the convergence of the pair 

The results of Theorem 2 may be paraphrased as follows: The mean- field distributions p^") converge weakly 
as measures on C(r2) (along a subsequence) to a Dirac mass concentrated at a minimizer of the Hamiltonian 
subject to the particle number constraint. In this limit, the variance of the field </>'■"■' vanishes while the variance 

in) 

of its gradient cpx is uniform over the interval SI, apart from the endpoints. This behavior is in good qualitative 
agreement with the long-time behavior observed in numerical simulations j^. 

Nonetheless, the continuum limit is degenerate in two ways. First, the limiting field (jf is deterministic. 
Second, the limiting gradient field, whose mean is 0* and covariance function is F*, is so rough that it does not 
even take values in the space of measurable functions (that is, it does not have a measurable version) [Q. These 
unusual properties of the limiting distribution are a direct consequence of the scaling of the inverse temperature 
/3'"' with n, which is required to maintain finite mean energy in the continuum limit. 



8 The Concentration Property 

We shall establish in this section the important result that the mean-field ensembles p*-"-*, which solve the 
maximum entropy principle (MEP), become equivalent in a certain sense to the microcanonical ensemble ( p^ 
in the continuum limit n oo. Specifically, we shall prove the following concentration property: 

Theorem 3. Let = 1, 2, • • •, be a sequence of solutions of the constrained variational principle (MEP). 

Then 

lim (iV„) = 7V° , lim Var(7V„) = , (56) 
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and 



lim (i7„> = H° , lim Var(iJ„) = . 



(57) 



This theorem has the interpretation that the mean-field ensembles p^") concentrate on the microcanonical 
constraint manifold Nn — N^,Hn — in the continumn limit. As we have emphasized above, it is a well- 
accepted axiom of statistical mechanics that the microcanonical ensemble constitutes the appropriate statistical 



equilibrium description for an isolated ergodic system pi] , 30 1. Thus, this concentration property provides an 
important theoretical justification for our mean-field theory. 

In the proof of Theorem 3, we will make use of the following elementary facts concerning Gaussian random 
variables: If is a Gaussian random variable with mean /i and variance (t^, then 



4r7 V . 



(58) 



We will also repeatedly make use of the result that l//?*^"^ = 0(n ^), which follows from equation (^ 



Proof of Theorem 3: The first conclusion in ( p6| ) is easy to establish. Indeed, using ( |2^ ) and the calculation 
(E^), we have 



k=l 



= N° + 0{n-^), asn 
To verify the second assertion in (|5^), we calculate using ( [l^ ) 



fc=i 



[Var(ufe) +Va,r{vk)] 



oo . 



Var(Ar„ 



\k=l , 
1 " 



(59) 



fc=i 



where, to obtain the second equality in this calculation, we have made use of the mutual statistical independence 
of the random variables mi, . . . , u„, wi, . . . , w„. But, as and Vk are Gaussian with means (ufc) and (wfe) and 
variances given by (0), we have from (|58| ) that 



(/3(")Afe)2 ^ /3(")Afc ' 
Upon substituting (p3) into (p9), we obtain 



hi)' 



4(«fc) 



(/3(")Afc)2 /3(»)Afc 



(60) 



Var(Ar„ 



«)^2 A? + ^(„) 



(/?( 

as n ^ cx) , 



Ai. 



since X]I-=i -^fc ^ converges as n ^ oo, and X]I-=i ^k ^((^fc)^ + ('^fc)^) — 2Aj ^A^° for all 
To establish the first part of (^7|), observe that from (p5|), there holds 

|(i?„)-i?°| = (e„(w("), «("))) -e„((u(")), («("))) 



(61) 



where the potential energy 0„ is defined by (17). Mimicking the calculations following equation (|26| ) and 
preceding equation (E9[), we find that 
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(e„(u("),z;(")))-e„((u(")),(«(")))| <C^[VarK)+VarK)] , 

k=l 

where C is a constant independent of n. Thus, according to (|6l|), (62) and (^5|), there holds 



(62) 



Next, we have that 



Var(ff„) = Var(i^„ + e„) < 2 Var(i\:„) + 2 Var(e„) , 



(63) 



where Kn is the kinetic energy defined as in (16). Using the statistical independence properties noted above, 
we calculate 



Var(i^„) = l/^(j2x,[(^ul-{ul)) + ivl-{vl))] 

k=l 

and using (|60|), we arrive at 

This last expression vanishes as n — > oo, because the sum X]fc=i ^k{{uk)'^ + (ffc)^), which represents twice the 
kinetic energy of the mean, is bounded independently of n. 

Finally, we show that the variance of the potential energy 0„ tends to in the continuum limit. Using the 
definition of 0„ and the Cauchy-Schwarz inequality, we have 




Var(e„) = - 



F((u("))2 + (t;("))2) _ (^((u("))2 + (t;("))2) 



dx 



< - Var (^^((m("))2 + (w("))2)) dx . 



(64) 



Now the potential F may be expanded about the mean ((u*^"^), (u'"))) as 



F((u("))2 + (i;("))2) = + (w("))2) 

+ 2/((w("V + (t.("))2) r(u("))(M(") - («("))) + («(") - («("))) 



(u(n) _ („("))) 



"^^^ >y („(") _ (^;("))) J 



where the matrix J is defined by ( p8| ) and (u("\ {;("') lies between (m("\d'^")) and ((u*-"^), (ij^"^)). Using the 
boundedness condition ^ on / and the statistical independence of u'-"^ and ti^"\ we obtain after a straightfor- 
ward analysis, that pointwise over U, 
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Var 



+ < C Var(ii(")) + Var(w(")) 



C 



- + - . (65) 



Here, C is a constant that is independent of both n and x € ^. As we have demonstrated in (p 



Var(«(")(x)) = Var(^(")(x)) = — ^ , (66) 

^ k=i 

and since u("^(a;) and u(")(a;) are Gaussian variables, it follows from ( |66|) and (|5^ ) that 



((^(")(a;) - {u(-\x))f) = ((z;(")(x) - (^("^(x)))^) = (^f] . (67) 

Thus, from (|6^)-(|67|), the fact that the eigenfunctions are uniformly bounded over independently of fc, and 
the fact that ^((il^))^ + (z;("))2) dx = 2iV" for all n, we have 

Var(6„) = 0(r7,^^) , as n ^ oo. 

This completes the proof of Theorem 3. 

9 Extension to Unbounded Nonlinearities 

In this section, we briefly indicate how our mean-field statistical theory, can be extended to a class of unbounded 
nonlinearities, which includes the focusing power law nonlinearities /(IV'P) = iV'l^'jO < s < 4. As illustrated in 
^, H, 1^] , numerical simulations for such power law nonlinearities exhibit the same phenomena as seen for bounded 
nonlinearities such as/(|V'P) = |'0P/(l + lV'P)j for which the results of numerical simulations are displayed 
in Figures (l)-(2). That is, the field ip approaches a long-time state consisting of a coherent soliton structure 
coupled with radiation or fluctuations of very small amplitude. As in the case of the bounded nonlinearities, the 
gradient ipx exhibits fluctuations of nonnegligible amplitude for all time. In fact, we expect, that this general 
behavior occurs as long as the nonlinearity / is such that the NLS equation (|l|) is nonintegrable and free of 
collapse (i.e., such that finite time singularity does not occur). Hence, we would like to apply our statistical 
theory to NLS with such nonlinearities, as well. 

Let us recall that to motivate the mean-field constraints (|2^ ) and (|2^) for nonlinearities / satisfying the 
boundedness condition (||), we invoked the vanishing of fluctuations hypothesis (|2ll). However, for nonlinear- 
ities such that the potential f (IV'P) grows more rapidly than Cj-i/'P as oo (as it does for the power 
law nonlinearities), the hypothesis (pi| ) is not sufficient to guarantee a priori that (0„(u'"', w*-"^)) converges 
)) as n — > cx). Thus the mean- field Hamiltonian constraint (|5|) can not be derived from 
the vanishing of fluctuations hypothesis (pT| ) alone for such /. Note that by making a stronger vanishing of 
fluctuations hypothesis, we could have weakened the assumptions on / and arrived at the same mean-fleld con- 
straints (|24|)-(p5|). In any case, we could simply impose these mean-field constraints and investigate the resulting 
maximum entropy ensembles p^"'\ If it can be shown that these ensembles exist and satisfy the concentration 
property expressed in Theorem 3, then we will consider this approach to be justified a posteriori. 

It is clear from Theorem 1 and Theorem 2 that, in order for our statistical theory to be well-defined, it is 
necessary that there exist minimizers (in iJp (f2), say) of the Hamiltonian H given the particle number constraint 
TV = iV°. This places restrictions on the class of nonlinearities that we can consider. However for nonlinearities 
/ such that the potential F satisfies 

There exists C > such that F{\i;\^) < C{\ip\^ + |V'|') , for some 2 < g < 6 , (68) 



2 



18 



it may be shown that such minimizers exist for any js^l- The condition (^Sj) is also crucial in establishing 
the well-posedness of the NLS equation (|^) as an initial value problem in the Sobolev space i?o(r2) Clearly, 
the focusing power law nonlinearities /(IV'P) = IV'I'* satisfy this condition when < s < 4. 

Thus, assuming that F satisfies the growth condition ( |68|), we may investigate the ensembles p^"^ which 
maximize entropy subject to the mean-field constraints (p^-(|25|). It is easy to see that the analysis of Section 
6 goes through without change. Thus, under the the maximum entropy ensemble /O*-"'', the Fourier coefhcients 
wi, • • • Un,vi, - ■ ■ Vn are mutually independent Gaussian variables, the mean-field minimizes the Hamiltonian 
given the particle number constraint 7V„ = A^", and the variances of the Uk and Vk are given by (^4|). It may be 
shown that the continuum limit of Section 7 and the important concentration property of Section 8 also hold. 



if we impose, in addition to (68), a lower bound on the potential F of the form 



There exists c > such that - c|V'|'' < F{\'ip\'^) , for some r > . (69) 



That is, Theorems 2 and 3 still hold for F satisfying ( |68| ) and (|69|), but the proofs are complicated by the 
unboundedness of /. Since our goal in this paper has been to emphasize conceptual issues rather than technical 
details, we shall not present the proofs here. However, we refer the interested reader to p7|], where some of the 



analysis has been carried out for nonlinearities satisfying (pq) and ( 69 ) 
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Appendix 

In this Appendix, we provide an analysis of the variational principle that defines the mean field corresponding 
to the maximum entropy ensemble p'"^. In particular, we prove that for (smooth) nonlinearities / satisfying 
the condition (|^), solutions of this variational principle exist, and we investigate the convergence properties of 
these mean fields as n ^ oo. 

First, we consider the variational principle 

7J(V') min , subject to e A° , (70) 

where the admissible set A° is defined as 

A° = {tA : -> C I V e H^{n) , N{tlj) = N"} . (71) 

The Hamiltonian H and the particle number N are defined by equations (||) and (^, respectively, and, as in 
the main text, Q — [0, L]. Recall also that the potential F is defined as F{a) — f{a') da' . 

Lemma Al. There exists a solution -0 g FIqI^I) of the variational problem ([70|). 

Proof: As / satisfies the boundedness condition (|^), there exists a constant K > such that l-FdV'P)! ^ -^IV'P- 
Thus, for all ^ G there holds 

HW = \ j \t.?dx-\ J Fm')dx 

> \j |^.|^dx~|7V0. (72) 
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Here, and in the remainder of this Appendix, all integrals are over the domain ft. It follows from (^2|) that 
H is bounded below on ^4°, so that we may choose a minimizing sequence {ipn)ne'N C A" for H. The in- 
equality ( [7^ ) also implies that the sequence -ipn is uniformly bounded in HQ{fl). Hence, by the Sobolev 
Imbedding Theorem [Q, there exists a subsequence, still denoted by ipm and a function ip such that ipn 
converges weakly to ip in -ffo(fi) and ipn converges in C{fl) to ip. Because of the convergence in C{fl), we 
have that / F{\ip\'^)dx — lim„^oo / i^d^/^np) dx. Also, since ipn converges weakly to i/j in HQ{n), there holds 
/ dx < liminfn^oo J \{ipn)x\'^ dx. Putting these two results together, we obtain that 

H{iP) < liminf = inf H{(l)) . (73) 



But, as '0n converges to '0 in (7(0), it follows that N{^) — N^, so that -0 G We conclude from (|73|), 
therefore, that ^ is a solution of the variational problem ( |70|) . 

Lemma A2. Any solution ip of the variational problem (^0|) is a solution of the ground state equation (^ for 
some real A. 

Proof: The conclusion of this lemma follows from an application of the Lagrange multiplier rule for constrained 
optimization. The parameter A in the ground state equation (H) is the Lagrange multiplier which enforces the 
constraint N{ip) = N°. 

We have demonstrated in Section M) that the mean fields ) corresponding to the 

maximum entropy ensembles p^"-* are solutions of the following variational problem 

71(0) min , subject to G A° , (74) 
where the admissible set A'^ is defined as 

AO ={0€S„|A^(0)=A^"} . (75) 

Here, S„ is the set of all complex fields on O having the form ip = u + iv, where the real fields u and v 
are elements of the n-dimensional space X„ spanned by the first n eigenfunctions ei, • • • , e„ of the operator 
—(P/dx^ on fl with homogeneous Dirichlet boundary conditions. S„ may be thought of as a closed subspace of 
the Sobolev space Hl{fl). The sequence of variational problems (Q) corresponds to a Ritz-Galerkin scheme for 
approximating solutions of the variational principle (|70|). The proof of the next lemma follows from arguments 
analogous to those that were used to prove Lemmas Al and A2. The conclusion about the smoothness of 
solutions of d?^ ) is obvious, because any such solution is necessarily a finite linear combination of the basis 
functions Ckix) = ■\/2/L sin(fc7rx / L). 

Lemma A3. For each n, there exist solutions of the variational problem (74). Any such solution 0^"-* is a 
solution of the differential equation 

0i^) + ^'^"^ (/(|0^"^ n0^"') + A(")0(") = , (76) 

where P*^"-' is the projection from -ffg (fi) onto S„, and A*^"-' is a Lagrange multiplier which enforces the constraint 
Ar(0(")) = N^. In addition, any solution 0;(") of (0) is in C°^{Vl). 

We wish to investigate the convergence of solutions of ( [74[ ) to those of (|7^). It is not our objective to 
obtain the strongest possible convergence results. For our purposes, it is sufficient to prove that any sequence 
of solutions of ( [z^ has a subsequence which converges in C^{fl) to a solution of ([70|). Such a result is all that 
is needed for the analysis of the continuum limit of the mean field ensembles in Section (|^). In general, we 
do not know that solutions of the variational principle are unique, so that subsequences can not be avoided. We 
shall now prove the following theorem. 

Theorem Al. If 0("), n = 1, 2, • • • is a sequence of solutions of the variational problem ([tJ) , then there exists 
a subsequence 0^"'^ of 0^^"^ which converges in C^[fl) as n' — > 00 to a solution of the variational problem (|70|). 
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Proof: Let -0^"'' be a sequence of solutions of and let F* = H{ip^"''>). Because A'^ c = 1,2, • • •, 

and A° = [J'^^i we have that H* > -ff*+i, and i?* \ H* as n ^ oo, where i?* is the minimum value of 
H over ^0. Thus, V^"^ is a minimizing sequence for the variational problem (|7^), and so following the proof of 
Lemma Al, there exists a subsequence -i/)^" ^ and a solution tp of the variational problem ( |70| ) such that i/;*^" ^ 
converges weakly in i?o(ri) to -0 and -0^" ^ converges in C(ri) to -0. 

Now, as -0*-" is a solution of the variational problem ([7^), it satisfies equation ( |76| ) for some Lagrange 
multiplier A'-" We shall now show that the sequence A*-" \n' oo is bounded independently of n'. For this 
purpose, we multiply equation (^6|) for -0*^" ^ by the complex conjugate -0'-"'^ and integrate over to obtain 

(77) 



It follows immediately from (^) and the fact that / |0>("')|2 rfx = 2iV" that 



|A("')| < 



2iV0 



p(«')(/(|^("')|2)^(«'))^(n')rfa; 



(78) 



Because ■0'^" ^ converges weakly to in i/o(Sl), the first integral on the right hand side of (fz^ ) is bounded 
independently of n' . The second integral on the right hand side of dTq) may be estimated as follows: 



p("')(J(|^("')|2)^("'))^(n')rfa; 



< 



< 



|0("')|2dx 
|0("')pdx 



|p("')(J(|^(n')|2)^(«'))|2^2. 

(/(|0^"'^l'))'lV^<"'^)l'rfa; 



< 27V«sup|/(|0("')n|. 

To obtain the first line of this display, we have used the Cauchy-Schwarz inequality. The second line follows 
from the fact that / |P'^")0p dx < J \(f>\'^ dx for all S H^{Q),7i e N, and to obtain the third line, we have once 
agam used the identity / |-0'"'^P dx = 2N°. Now, owing to (|^), sup2.gj^ 1/(10''^" is bounded independently 
of n' , so the preceding calculation demonstrates that the second integral on the right hand side of ( |78| ) can also 
be bounded independently of n' . This implies that A*^" ' is uniformly bounded in n' . 

The strategy now is to use the uniform boundedness of the eigenvalues A^" ^ to establish that the subsequence 
^ is uniformly bounded in the Sobolev space _ffg(ri). From this result and the Sobolev Imbedding Theorem 
iQ, we may conclude that the 0''-" converges in C^{n) to the function -0 as above, which is a solution of the 
variational problem (kO|). This is the desired conclusion. 

To prove that 0;^^ is, in fact, uniformly bounded in Hg{fl), we multiply the equation ( [76| ) for 0;*^" ^ by 
the complex conjugate of (0^" ^)xx and integrate over Q. This yields, after integrating by parts and using the 
homogeneous Dirichlet boundary conditions, 

|(V'("'^)..pdx = A("') / |(^("')),|2dx+ /^,(p("')(/(|V'("'f )0'("'))).rfa;. 



(79) 

Using the Cauchy-Schwarz inequality, and the inequality J \{P"'{(j>))x\'^ dx < /|0a;pc?a;, which is valid for all 
n G N and all (j) e Hq (Q) , we find that 



0;(«')^(p(«')(/(|V,("')|2)0("'))),da; 



< 



' dx 



("')|2A,/,("')V.|2 



i(/(i^^"r)V' 



dx 



But, a straightforward calculation yields that, pointwise on fl, 

i(/(iV'("')nv;("')),|2 < 2(/(iv.("')n)2|(^(«')),|2+8(/'(iV'("'f ))2|^("')n(0'("')).i^ 
< q(^^"'^).p, 



(80) 



(81) 
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for a constant C independent of n' . The second line of this calculation follows from the boundedness condition 
(|) on /. Taken together, @-(|l|) imply that 

J < (C+ |A("')|) J dx . (82) 

As we have established above, A'" ^ and / KV'*-" dx are bounded independently of n' , so that from (^2|), we 
conclude that the subsequence V'^" ^ is uniformly bounded in ifo(^)- This concludes the proof of Theorem Al. 
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Figure 1: Evolution of \ip{x,t)\'^ for tlie NLS equation iifjt + ijjxx + V'|V'P/(1 + l"^!^) = 
on the spatial interval [0,64], with homogeneous Dirichlet boundary conditions imposed. The 
initial condition is tp{x,0) = 0.5 tanh(x) tanh(64 — x) plus a small random perturbation. The 
integration scheme used is the split-step Fourier method with 512 Fourier modes. The relative 
deviations {N{t)—N^) /N^ and {H{t)—H^)/H^ of the particle number N{t) and the Hamiltonian 
H{t) from their initial values and are less than 3 x 10~^ and 8 x 10~^, respectively, for 
the duration of the simulation, (a) t=100; (b) t=1100; (c) t=100100; (d) t=200100. 

Figure 2: Evolution of \ipx{Xit)\^ for the same equation, boundary conditions, and initial 
condition as in Figure 1. (a) t=100; (b) t=1100; (c) t=100100; (d) t=200100. 
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